A systematic study of direct photon production in heavy ion colhsions 
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A theoretical derivation of photon bremsstrahlung, induced by the interactions of an energetic 
quark in a hot and dense quark-gluon plasma, is given in the framework of the reaction operator 
approach. For the physically relevant case of hard jet production, followed by few in-medium inter- 
actions, we find that the Landau-Pomeranchuk-Migdal suppression of the bremsstrahlung photon 
intensity is much stronger than in the previously discussed limit of on-shell quarks and a large 
number of soft scatterings. This result is incorporated in the first systematic study of direct photon 
production in minimum bias d-|-Cu and d-|-Au and central Cu-|-Cu and Au-|-Au heavy ion collisions 
at the Relativistic Heavy Ion Collider at center of mass energies ^/s = 62.4 GeV and 200 GeV. We 
find that the contribution of the photons created via final-state interactions is limited to 35% for 
2 GeV< pt < 5 GeV and at high transverse momenta the modification of the direct photon cross 
section is dominated by initial-state cold nuclear matter effects. 

PACS numbers: 12.38.Bx,12.38.Mh,24.85.+p,25.75.-q,25.75.Cj 



I. INTRODUCTION 

In the highly successful hard probes program at the 
Relativistic Heavy Ion Collider (RHIC), the interplay of 
nuclear effects that alter the cross section for direct pho- 
ton production is not yet well understood. Direct 7 mea- 
surements [l|, 0] have provided an important baseline to 
help establish the dominance of final-state effects for the 
observed large suppression of energetic hadrons [1, 0, [H| 
in nucleus-nucleus collisions. The interest of the theo- 
retical community, however, has been centered on possi- 
bly large new sources of photons as a by-product of an 
energeticparton propagating in the quark-gluon plasma 
(QGP) Early work predicted significant enhance- 

ment of low to intermediate px photons 0, Q based 
on a 7 emission pattern characteristic of on-shell quarks 
and a very large number of soft interactions in matter, 
e.g. 0]. This limit yields only a small, ^ 30%, Landau- 
Pomeranchuk-Migdal (LPM) suppression of the photon 
spectrum relative to the incoherent Bethe-Heitler result. 
It does not reflect, however, the possibly large cancella- 
tion between the bremsstrahlung associated with hard jet 
production and the subsequent quark scattering, which 
was found to control induced gluon emission in finite non- 
Abelian plasmas [1, [13, [HI ■ Inverse Compton scattering 
in the QGP [l^] has also been suggested as the domi- 
nant source of intermediate pt photons if the 7 acquires 
most of the momentum of the incoming quark. Such en- 
hancement effects, in conjunction with their associated 
negative azimuthal asymmetry coefficient V2{pt) fsl. [13|. 
do not appear to be compatible with existing direct pho- 
ton data at RHIC [2,[i3,[ll- 

Recent phenomenological refinements [IB, [i3l suggest 
that the QGP enhancement of direct 7 production may 
be smaller than previously expected and partly can- 
celled by the quenching of fragmentation photons. Still, 
there is no calculation to date that consistently includes 



known nuclear matter effects, such as the Cronin ef- 
fect 0, [H, [13, shadowing [Hi, [12, and cold nuclear 
matter energy loss |23j , to provide quantitative guidance 
for the relative strength of initial- and final-state modi- 
fications in the observed 7 cross sections. An additional 
serious deficiency in the theory and phenomenology of 
direct photon production in heavy ion collisions is the 
absence of systematic studies in proton-nucleus (p-|-A) 
and nucleus-nucleus (A-l-A) collisions for different system 
sizes and center of mass energies. This is especially true 
now, when new experimental results from RHIC are soon 
expected to become available [13, [IS]- Last but not least, 
only through extensive detailed comparison between the- 
ory and data [^ can one test the model validity and gain 
confidence in the extracted quantitative properties [25| 
of the dense matter created in heavy ion reactions. 

With this motivation, we derive the QGP-induced 7 
spectrum for hard quark production in finite size plas- 
mas. The same model of jet-medium interactions is used 
to calculate the quark conversion cross section and the 
suppression of fragmentation photons. These theoreti- 
cal results, when applicable, are combined in a numer- 
ical simulation with cold nuclear matter effects to pro- 
vide model predictions for d-|-Cu, d-|~Au, Cu4-Cu and 
Au-|-Au reactions at center of mass energies of 62.4 GeV 
and 200 GeV per nucleon pair at RHIC. This article is 
organized as follows: in section [IT] we highlight the dif- 
ferences between gluon and photon bremsstrahlung and 
identify the theoretical approach that reproduces the 
known Bethe-Heitler spectrum. The derivation of the 
final-state medium-induced photon radiation in the reac- 
tion operator approach is given in section [Till Numerical 
results, relevant to the phenomenology of direct 7 pro- 
duction are also shown. In section IIVI we carry out a 
systematic investigation of cold and hot nuclear matter 
effects that alter the mid-rapidity photon cross section 
in ultra-relativistic collisions of heavy nuclei at RHIC. A 
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summary and conclusions are presented in section IVl 



II. PHOTON VERSUS GLUON 
BREMSSTRAHLUNG 

The computation of photon bremsstrahlung is usuaUy 
considered to be easier than that of gluon bremsstrahlung 
due to the absence of self-interactions of the gauge boson. 
However, it is not well appreciated that the two physics 
processes are quite different. To illustrate this, we first 
examine the radiative amplitude that corresponds to the 
case of single scattering of a fast on-shell quark [IB] : 
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In Eq. ([!]) gs is the strong couphng constant, = 
[fc"*", k_L] is the momentum of the radiated gluon in 
light-cone coordinates and e is its polarization vector. We 
denote by = [q~^ = 0, q~ ,q±] the momentum exchange 
with the medium at position z and by T^, T" G SU(3) the 
color matrices at the emission and interaction vertexes. 
Evidently, it is the color rotation of the parent parton and 
the re-interaction of the bremsstrahlung gluon in nuclear 
matter that determine the gluon emission intensity and 
allow neglect of the deflection of the jet. If we, how- 
ever, take the Quantum Electro-Dynamics (QED) limit 
gs ^ eq ^ (±l/3,±2/3)e, ^ 1 in Eq. Q we find 
M.1^^{k) 0. Therefore, a theoretical approach devel- 
oped to describe gluon emission cannot be directly gen- 
eralized to photon emission and vice versa p7| . Our con- 
clusion is independent of the specific example of incoher- 
ent parton scattering. All regimes of coherent inelastic 
scattering can be treated in the unified framework of the 
reaction operator approach [23j . The reaction operator 
Rn describes the effect of one additional correlated in- 
mcdium scattering at position z„ at the cross section level 
and is process dependent. Taking the above mentioned 
QED limit of Rn derived for gluon bremsstrahlung [T^], 
one finds: 



Rn 
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In Eq. ^ Cp is the quadratic Casimir in the fundamental 
representation of SU(3). Thus, better treatment of jet- 
medium interactions is needed for both incoherent and 
coherent photon emission calculations. 

With these results in mind, we first identify the re- 
finement of the kinematic approximations necessary to 
derive the induced 7 spectrum. The scattering of a fast 
quark in nuclear matter is modelled via interactions with 
an external non-Abehan field V^-''^{q) [ist : 



V^^'{q) = n^27rS{q+)V^{q)e"'-' 
gsV^q) = v{q)T^t) . 



(3) 



Here, the four- vector = 6^^' = [0, 1, Oj^] and the color 
matrix T^(t) £ SUc(3) in Eq.([3|) represents the target 



charge that creates the non-Abelian field. We take the 
Fourier transform v(q) to be of color-screened Yukawa 
type but with Lorentz boost invariance: 

_ Anas Anas , ^ , 

where we have used the (7+ = choice of frame. This 
specific form of v(ci) = v{(\i_) is particularly useful since 
in-medium interactions in both hot and cold nuclear mat- 
ter are of finite range Tint. = and we shall assume 
that Ag/z 3> 1, where A is the quark mean free path. 

The differential photon bremsstrahlung spectrum 
arises from single-Born scattering diagrams shown in the 
top panel of Fig. [1] Using a high energy approximation 
for the quark to simplify the interaction and emission 
vertexes we obtain: 
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For the second 7 emission diagram similar considerations 
lead to: 
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and the radiative matrix element at position Zi reads: 



Mradik,{i}) = e 



e-pf €-pi 



k ■ Pf k ■ Pi 
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(7) 




V(qi). V(q)) 



V(qi). V(q)) 



FIG. 1: Top panel: single-Born diagrams formedium-induced 
7 emission. Bottom panel: the corresponding double-Born 
diagrams in the Zii — > Zi limit. The third diagram, known to 
vanish in this limit [l3| , is not shown. We denote by " x " the 
propagators that enter the calculation. 
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In Eq. ([7]) the coUisional amplitude is not shown. 
Let the initial- and final-state momenta of a fast on- 
shell quark be ^ [E+ , Ql ^_J{2E+)M^ t-i], Pf = 
[E+, Ql J{2E+), such that , - = ,. 

The double differential medium-induced photon distribu- 
tion is then given by: 



,^^dN'<{k;{{\) 
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and is dominated by emission coincident with the di- 
rections of the incoming and the outgoing quarks. 
Changing variables k = kx — k|j'|'°''''' , k^"'''^ = 
Q± ik^ / E'^ ,Q± i_ik^ / E'^ , respectively, we obtain the 
QED double logarithmic result: 
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k+ 



■In- 



(9) 



where regulates the coUinear divergence. 

In the case of coherent gluon emission in finite me- 
dia with few subsequent scatterings, the interference be- 
tween the hard (vacuum) and soft (medium-induced) 
bremsstrahlung largely determines the LPM cancellation 
pattern 0, [lol |. Double-Born diagrams, with two mo- 
mentum exchanges at the same position Zi ~ z[, can 
contribute at any fixed order in opacity. The relevant 
double-Born diagrams for photon emission are shown in 
the bottom panel of Fig. [1] We find: 
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To obtain the result in Eq. pOj) we averaged over the 
initial and summed over the final parton colors. We also 
carried out the average over the position of the scatter- 
ing center in the transverse plane: J d'^'z± exp[— izj^ • 
(q-L + ql)] - Al\2TT)^S^{qi_ + q'^). For single-Born 
interactions, see Eqs. ([5]) and (O, such averages are 
possible only after squaring the amplitudes. The result 
differs from the one for virtual interactions in yielding 
(5^(qx — q'j^), and in the absence of the factor —1/2, which 
accompanies the collision term. 

Similarly, for the second virtual diagram in Fig. [1] we 



obtain: 



A, 



X e — 



1 rfqj_ 

2 J (27r)2 
e-p, 



Hq^)\ 



k -pi 



(11) 



Adding Eqs. (fTU]) and (fTTI) we see that the same gen- 
eral radiation matrix element, Eq. can be factor- 
ized for double-Born interactions. In Eqs. (fTO]) and (fTTj) 
dcr°Vd2q_L = (CtCp/dA)|w(q±|V(27r)2 is the differential 
scattering cross section, calculated in the Born approxi- 
mation. However, qj_-l-q'j_ = 0, andp,; = Pf implies that 
the double-Born interaction does not contribute a new 
photon bremsstrahlung amplitude. 



III. DIFFERENTIAL PHOTON 
BREMSSTRAHLUNG SPECTRUM TO ALL 
ORDERS IN OPACITY 

To define the iterative procedure of computing the pho- 
ton bremsstrahlung contribution from multiple scatter- 
ing, we first consider the action of the direct operator 
Dn at position z„ on a radiative amplitude with n — 1 
correlated scatterings. In what follows we have dropped 
the collisional amplitudes since they were shown to yield 
an elastic scattering cross section per order in opacity for 
both diffusion [H, [2^ and radiative [l^l processes. The 
result of such action can be represented as: 

A.x™'.„_,(fc) ^ (i + s„)A<™.^.,_^(fc) 

= ■M™.'..„_,(fc) + (-i )'^"^^'--'"-^ A^™d(fc, {n}) . (12) 



Here, the factor 



arises because every virtual con- 
tact interaction in the amplitude gives a factor — i, and 
Ny{Ml^^i^_J is their number. The first term in Eq. ^ 
corresponds to a momentum exchange with the energetic 
jet. In the high energy limit we do not keep track of 
the transverse modification of the parent parton except 
for the contribution to the soft photon bremsstrahlung. 
This approximation does not affect the calculation of the 
intensity spectrum dl^ /dk'^ but the angular distribu- 
tion must be convoluted with the medium-induced jet 
acoplanarity, which we here neglect. The second term 
in Eq. (I12p is the new radiative contribution, Eq. ([7]), at 
position z„, and the prefactor accounts for the number 
of preceding virtual interactions in the amplitude. Next, 
we consider the double-Born interaction of the quark at 
position z„. From section|lI]wc know that there is no new 
7 bremsstrahlung contribution since pf = Pi- However, a 
factor —1/2 arises for the forward elastic scattering. The 
modification of the amplitude _^{k) is found to 

be: 
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In both Eq. (|12p and Eq. ([13]) wc have omitted the color 
factors since for the simple case of individual parton prop- 
agation these are trivially absorbed in the elastic scatter- 
ing cross section or inverse mean free path per order in 
opacity (isj . 

Let z„ = 0, 1, 2 indicate for each n: no interaction, di- 
rect interaction and virtual interaction with the medium, 
respectively [loj . Expansion in powers of cr"' or, cquiv- 
alcntly. 1/Xq requires that in the conjugate amplitude 

-^rad*"(^) ^® have in = 2 — in- Consequently, the con- 
tribution to the radiation pattern at n-th order in opacity 
is: 



"ii-'-^n-i— 
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= E -^taT^""^ {k){BiBn + Bi + S„)A^™.^.._, (fc) . 

il-'-in-l—O 
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For medium induced photon emission, the reaction op- 
erator i?„ = DlDn + V,l + Vn= BlBn + B,t + Bn has 
a particularly simple form. The first term in Eq. (jl4|) 
vanishes beyond first order [n = 1) in opacity [loj : 



where we have used (— ^ — 5 + l)" ^ = for 71 > 2 p^ . 
The 2"^^ and Z''^ terms in Eq. (Ull) yield: 



theoretical result derived in this letter reads: 
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if ri > 3. Therefore, unlike the case of giuon 
bremsstrahlung, photon bremsstrahlung contributions 
vanish beyond second order in opacity. 

Taking into account the interactions of the parent 
quark along its trajectory through the QGP and the av- 
erage over the transverse momentum transfers, the main 
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In Eq. P?)l Azf = zf — zf_-^^, the Az; integrals are nested, 
and tJ'^ = \s?/{2u) « ^2 k" is the inverse photon forma- 
tion time. When Tj^Xq 3> 1 the photons decohere early 
from the parent quark and our result reduces to inco- 
herent emission from individual scattering centers. Our 
result differs from previous findings 0, 0] in several im- 
portant ways: first, it treats the case of finite L/\q ^ few, 
relevant to heavy ion physics. Second, we find that the 
most significant contribution to the LPM effect for pho- 
tons comes from the interference of the medium-induced 
photon radiation with the hard emission from the large 
scattering of the parent quark. Finally, we find that 
there can be non-linear corrections ~ to the domi- 
nant linear in L behavior of the photon spectrum. Note 
that calculations of final-state 7 emission have also been 
carried out in deep inelastic scattering on nuclei [13, 
using the high twist approach [ll[. 

We are now ready to study numerically the final-state 
QGP-induced photon spectrum. Our results arc limited 
to first order in opacity. The question of whether de- 
structive interference may even cancel part of the hard 
photon bremsstrahlung is deferred for future studies. 
With Eq. ([7]) representing the amplitude of both light 
and heavy fermions, the intensity spectrum is easily 
generalized to quarks with physical and thermal mass, 
771^ = A'/2 ^ CpglT^/^:. Notable differences from the 
TO = case include regulation of the poles in Eq. ([5]) via 
(k'^ / E'^)'^m? , appearance of new terms oc m?, and finite 
quark velocity, < 1. Wc show an example of a quark 
propagating outwards from the center of the medium cre- 
ated in 6 = 3 fni Au+Au collisions at RHIC. For a typical 
integrated gluon rapidity density, dN^/dy ~ 1150, dis- 
tributed proportional to the 2D participant number den- 
sity diVpart /rf^x^ , we calculate the temperature T{x±,t) 
clS cL function of time and position in the transverse 
plane, assuming approximate Bjorkcn expansion of the 
QGP near midrapidity y = 0. The necessary Debye 
screening scale ^J.{T) ~ m]j(T) = gsT, elastic scat- 
tering cross section cr'3(T), and quark mean free path 
A,(r) = l/c7i9(T)p{T) are evaluated as in [3 with = 
2.5. The top panel of Fig. [5] shows the medium-induced 
photon number spectrum dN'</dx = {e/ Cq^dN^ / dx, 
X = k^ /E'^, normalized by the squared fractional quark 
electric charge (bar will denote such scaling for any 
physics quantity). We considered light, M^^d = GeV, 
and heavy, Mc = 1.5 GeV, Mi, = 4.5 GeV, quarks of en- 
ergy Eq ~ 100 GeV in the Bethe-Heitlcr limit and to first 
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FIG. 2: Top panel: medium-induced photon number spec- 
trum versus x — jE'^ for Eq = 100 GeV light, charm and 
bottom quarks in central Au-|-Au collisions at = 200 GeV. 
Both the incoherent Bethe-Heitler (solid) and the coherent 
final-state (dashed) bremsstrahlung cases are shown. Bot- 
tom panel: partial fractional quark energy loss AiSg^'/Sq 
due to 7 emission versus the mean number of quark inter- 
actions (ji) = L/\q. Insert illustrates the energy dependence 
of AEi~'^/Eq. 



from the naive expectation, AEq oc Eg, as illustrated 
in the insert in Fig. [21 Medium-induced photon emis- 
sion in both the coherent and incoherent limits is too 
small, I\E^'^ /Eq < 1%, to contribute to the quenching 
of quark jets. Last but not least, we emphasize that the 
number of interactions in the QGP, even for jets emerg- 
ing from the center of the heavy ion collision region, is 
small, {n) = 3.1. As in the case of gluon emission, this is 
a clear indication that, even in the most central Au+Au 
reactions at RHIC, we are not in the limit of large num- 
ber of scatterings. For consistency with the calculation 
of light hadron attenuation, in our numerical estimates 
we use an effective geometry of L = 6 fm and uniform 
distribution of partons [J. In this case (n) is reduced 
to 2.4. accounting for jets close to the periphery of the 
interaction region. 



IV. PHENOMENOLOGY OF HARD PHOTON 
PRODUCTION IN P+A AND A+A COLLISIONS 

With the results from the previous section at hand, 
we now turn to the question of hard, pT > a few GeV, 
photon production in heavy-ion collisions. QGP and cold 
nuclear matter effects are identified through the nuclear 
modification ratio: 

where the number of binary collisions. N'^g{b) in 
Eq. p^ . is computed in an optical Glauber model. The 
baseline p+p cross section is evaluated in factorized per- 
turbative QCD to lowest order and leading twist as fol- 



order in opacity. In the absence of coherence, the induced 
7 spectrum is dominated by photon energies uj ^ a. few 
xT. The contributions of the heavy quark sector to the 
medium-induced photon multiplicity and the energy loss 
due to photon emission are strongly suppressed. When 
the interference between the vacuum and the medium- 
induced photon radiation is taken into account, we find 
that the spectrum dN^ /dx is suppressed much more ef- 
fectively than in the limit of very large number of soft 
scatterings for on-shell quarks 0]. The bottom panel of 
Fig. [2] shows the dependence of the partial fractional en- 
ergy loss due to 7 emission versus the mean number of 
quark interactions in the QGP: 

(n) = L/\ = J dz/\{z) . (18) 

We note that for light and moderately heavy quarks, such 
as the charm quark, there is clear non-linear dependence 
of AE^''^/Eq on (n). It is correlated with deviations 



10-' 




Pt [GeV] 



FIG. 3: Direct photon production cross section in p-|-p col- 
lisions at ^ = 62.4 GeV and 200 GeV. Data at the higher 
RHIC energy is from PHENIX [sj . Insert shows the fraction 
of fragmentation to all direct photons. 



6 



lows: 



fih)f{h) 
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-^D^/c{z,fif 



2S 

, 4>a/N {Xg , Hf ) (/)b/Af {Xb , fJ.f ) 
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Mab^cd\ 



(20) 



In Eq. PU| we adhere to the standard notation, see for 
example P. [2ll. [33| . and the Jacobian reads: 



s 

,xtika, kb) = — 



1 



(21) 



subject to the hard scattering constraint ka^b < Xa,b^fS- 
In elementary p+p collisions we use (fc^ = 0.9 GeV^ for 
the normalized Gaussian parton transverse momentum 
distributions, f{ka.b), and phenomenological K-factors 
cancel in Rab{pt)- One observes that Eq. (^0)) includes 
both prompt photons, c = 7, ac = (Xem = 1/137.036, 
D^/^[z) = (5(1 — z), and fragmentation photons, c = 

M~, (7, ac = as(/ir), with D^/^{z, ^fr) taken from Ref. 
. Figure [3] shows the cross sections for direct 7 pro- 
duction in p+p coUisions at ^/s = 62.4 GeV and 200 GeV 
and identifies the fragmentation fraction of direct pho- 
tons that will undergo final-state modification in the 
QGP. 

The first and necessary step, that has thus far been ne- 
glected in direct photon phenomenology, is a systematic 
study of nuclear effects in p-l-A reactions. Not only do 
the isospin effect, Cronin effect, shadowing, and energy 
loss in nuclei significantly affect the observable Raa{pt), 
but these are interesting in their own right in light of the 
current and future d-|-Au measurements at RHIC [2^ . 
Theoretical approaches to the Cronin effect are well doc- 
umented [1^ and are centered around the idea of initial- 
state multiple scattering [23|. In our calculations the 
Cronin effect in the nuclear medium is controlled by the 
broadening of the transverse momenta of incoming par- 
tons as follows: 



{k 



t) fned 



X 



X max[l,ln(l + (5p|)] , (22) 



g,9 



where the term ln(l -I- Sp^) takes into account the power 
law tail of the k-r distribution beyond the naive Gaus- 
sian random walk [33| (e.g. Moliere multiple scattering). 
Specifically, = o.l2 GeV^, Ag = (Cj./Ca)A, = 1 fm 
and S = 0.3 arc fixed to reproduce the experimental p+A 
data when wc include initial-state energy loss and arc 
compatible with our earlier findings for {k'^)med within 



20% [20|, l30|. In Eq. ^ L gives the length of the nu- 
clear medium. We note that in semi-inclusive deep in- 
elastic scattering (SIDIS) the final-state broadening is 
expected to be approximately equal to the initial-state 
diffusion in fcy. A naive estimate of the medium-induced 



contribution to the hadron {Ap'^)a = {Pt)a ~ {Pt)d 
can be made from Eq. (HH). The numerical results at 
leading order for ^ 5 GcV^, {Ap^)Ne = 0.14 GcV^ 
{Ap^)Kr = 0.29 GeV2 and {Ap^)xe = 0.36 GeV^, 
are compatible within error bars with recent preliminary 
HERMES data HH. 

In this work, only the EMC effect was 
parametrized (2^ and shadowing was calculated from 
the coherent final-state parton interactions [2l|. The 
only parameter, the scale of power corrections [2l|, is 
constrained by the mean squared momentum transfer per 
unit length /i^/A so that (C^A^/^)^ ^ ~ [2y?L/X)q^g in 
minimum bias reactions, which yiel ds ( £^)o w 0.12 GeV^. 
Cold nuclear matter energy loss [23| was consistently 
calculated with the same lengths, momentum transfers 
and mean free paths as in Eq. (|22p and incorporated in 
the pQCD calculation as follows: 

'l^a.b/N{Xa.,b,Q'^) (l)a,b/N ^ ^ ^^'^ ^ , Q^^ 

~ (l)a,b/N {xa,b{^+'^a-b),Q'^) , (23) 

when ea,b 1- Here, f, arc the fractional energy 
losses for the incoming partons a, b evaluated in the 
rest frame of the corresponding target nucleus. It has 
been shown for the typical hot QGP created at RHIC 
that the effect of multi-gluon fluctuations of the induced 
bremsstrahlung on the quenching of jets could be mim- 
icked by a simple reduction of the mean energy loss 
Ceff = kAE/E, where k = 0.4 - 0.5 Though fluc- 
tuations of the initial-state cold nuclear matter energy 
loss have not yet been studied, in our calculations we 
make a similar approximation. We note that the weaker 
falloff with energy of the incoming parton flux, when 
compared to the hard scattered (additional ^ ^/Pt) final- 
state quark and gluon distributions, implies values of k 
closer to unity. The specific choice that we make, k = 0.7, 
was constrained through the evaluation of the tt^ produc- 
tion cross section in Cu+Cu collisions at RHIC Q which 
provided a good description of the experimental data. 

Our results for minimum-bias d-f-Au (solid lines) and 
d-|-Cu (dashed lines) collisions are shown in Fig. U) The 
calculated RdAipr) at = 200 GeV and 62.4 GeV are 
presented in the top and bottom panels, respectively. In 
the pt < 5 — 7 GeV region, the nuclear modification 
is dominated by Cronin enhancement with a magnitude 
that is not strongly affected by shadowing but is sensi- 
tive to the initial-state energy loss. The complementary 
PT > 5 — 7 GeV part of phase space is characterized 
by RdAiPr) < 1 with the isospin effect, included in all 
calculations, being one of the major contributors. The 
EMC effect only becomes noticeable at the lower CM. 
energy and at the highest transverse momenta. Initial- 
state, cold nuclear matter energy loss can contribute as 
much as the isospin effect at high px- We emphasize that, 
at transverse momenta 15 GeV, nuclear effects on the 
direct photon cross section can be as large as 20% at 
= 200 GeV and 40% at = 62.4 GeV in minimum 
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FIG. 4: Cold nuclear matter effects manifest in direct pho- 
ton RdAuijPr) in minimum bias d+Au (solid lines) and d+Cu 
(dashed lines) collisions. RHIC energies ^/s — 62.4 GeV and 
200 GeV are shown in the top and bottom panels, respec- 
tively. Preliminary — 200 GeV minimum bias d+Au data 
is also included [2411. 



bias d+A collisions. Preliminary experimental data 24 1 
is consistent with the theoretical expectation. However, 
given the large error bars, it cannot discriminate be- 
tween different cold nuclear matter effects or constrain 
their magnitudes. Careful experimental investigation is 
needed to pinpoint these effects and, since they will be 
enhanced in A-|-A collisions, caution should be exercised 
in the interpretation of the R\j^{pt) findings. 

In A-l-A collisions, QGP-induccd modification of di- 
rect photon production cross section includes compet- 
ing effects: the quenching of fragmentation photons and 
the tree level (jet conversion) and bremsstrahlung pho- 
ton enhancement. The quenching of fragmentation pho- 
tons is modelled in the same way as the quenching of 
hadrons [20| and can be combined with the medium- 
induced 7 contributions as follows: 



l-z 



de P{e) 
dz 



1 



-D. 



+ ^elnv.(c)'5(l-z). (24) 



Here, P(e) is the probability distribution of the fractional 
jet energy loss e = I^E/E 0, HI], dN^^^ {c)/dz is the 
QGP-induced bremsstrahlung that we calculated in sec- 
tion |TT1 and N^onv (c) is the number of jets converted to 
photons. For such conversions, in Eq. (j24p the py « pc 



approximation [l^l that stems from the limit of small t- 
channel momentum transfers to energetic jets is implicit. 
For a parent quark propagating through the QGP: 



ATT 

conv, 



(c) 



dt piT)aZ- 



■79 



(T) 



(25) 



where time and position dependence is taken via 
r(xx,t). The cross section in Eq. (pS)) . with s ~ 2mDE 
and t e (to£), s/4), consistent with the forward scattering 
approximation, reads: 



■nasaem , E 
■ In ■ 



(26) 



Numerical results in central Au+Au and Cu+Cu col- 
lisions at the intermediate and top RHIC energies are 
shown in Fig. |5l The final-state gluon rapidity densi- 
ties dN9/dy = 1150 {Au, 200 GeV), 800 {Au, 62 A GeV), 
370 (Cu,200 GeV), 255 (Cu,62.4 GeV) are constrained 
by measured and extrapolated particle multiplicities at 
RHIC [|. Of the final-state effects in the QGP, the left 
panels only include the quenching of fragmentation pho- 
tons. Since a small fraction of the direct 7 come from 
fragmentation processes and quark attenuation is signif- 
icantly smaller, ~ Cf/Ca, when compared to gluon at- 
tenuation, the observable QGP modification is also very 
small. In fact, Raa{pt) is dominated by cold nuclear 
matter effects, such as the ones shown in Fig. 01 ampli- 
fied by thepresence of two large nuclei. Raa{pt) from 
published [J, [s^l and preliminary 0, [l3l data are shown 
for comparison. Given the error bars, only large Cronin 
enhancement at y/s = 62 GeV is excluded, suggestive of 
the role of initial-state inelastic jet scattering in control- 
ling the magnitude of the Cronin effect. 

The right panels in Fig. |5| include the medium- induced 
photons and the jet conversion contribution. For com- 
pleteness, we have also shown results for the incoherent 
Bethe-Hcitler radiation spectrum and absence of initial 
state energy loss. This scenario gives direct 7 enhance- 
ment as large as factors of 3 and 4 over the p+p baseline 
at = 200 GeV and 62 GeV, respectively. Not only 
is this case not supported theoretically by the results 
derived in this paper and in Ref. psj . but when com- 
pared to data, even with the present large experimental 
error bars, it is clearly excluded. In calculating the co- 
herent final state photon emission rate, Eq. pT]) . and 
the jet conversion rate, Eq. (|25|) . we also account for the 
time-dependent quenching of the quark as it propagates 
through the QGP: 



Raa{pt, t) = {l- f{t)) + RAA{PT)f{t) . 



(27) 



Here, Raa{pt) is the full final-state quark quenching, 
evaluated as in Q for different system sizes and center 
of mass energies, and f{t) interpolates between and 1 
to give the time dependence of radiative energy loss [2^ . 
We find that this effect reduces dN^^Jdz and N^o^v. t>y 
30%. Our results show that at transverse momenta 
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FIG. 5: Calculations of direct photon Raa{pt) in central Au+Au (solid lines) and Cu+Cu (dashed lines) collisions. RHIC 
energies ^/s = 62.4 GeV and 200 GeV are shown in the top and bottom panels. Of the final-state QGP effects left panels only 
include the quenching of fragmentation photons. Right panels also include the jet conversion and medium-induced y ( both 
coherent and incoherent limits) enhancements. Preliminary = 200 GeV and 62.4 GeV Au+Au data is shown 0, jlj- At 
200 GeV two preliminary Raa{pt) results from P, [sj (circles) and (squares) are extracted. 



Pt < 5 GcV the contribution from jet conversion to the 
total photon cross section is limited to ~ 25% and the 
contribution of medium-induced 7 is limited to ^ 10%. 
In the high pr range the total enhancement contribution 
is found to be ^ 5%. 



V. CONCLUSIONS 

In this paper, we provided a theoretical derivation 
of the final-state QGP-induced photon bremsstrahlung 
for the experimentally relevant case of hard jet produc- 
tion. We demonstrated that while the physics processes 
that control photon and gluon emission differ, the com- 
mon Landau-Pomeranchuk-Migdal interference between 
the radiation from the hard scattering and the radiation 
from the subsequent soft quark interactions in the plasma 
leads to a significant suppression of the 7 intensity. The 
photon spectrum was found to be attenuated at least by 
a factor of several for jet energies relevant for RHIC phe- 
nomenology, in contrast with the estimated modest 30% 
attenuation for asymptotic t = —00 on-shell jets in the 
limit of a very large number of soft interactions, where 
the interference with the bremsstrahlung that accompa- 
nies hard jet production has not been taken into account. 
We found that the suppression of /dk~^ also leads to 
non-linear dependence of AE'^{L/Xq,Eq) on the system 
size and sub-linear dependence on the parent quark en- 



ergy that have previously been associated primarily with 
gluon emission. 

To help identify the significance of both cold and hot 
nuclear matter effects on direct photon production we 
carried out the first systematic phenomenological study 
of R\b^Pt) in midrapidity d+-Cu, d+Au, Cu+-Cu and 
Au-|-Au reactions at RHIC energies of ^/s ~ 62.4 GeV 
and 200 GeV. As expected, in all cases we found that in 
the absence of QGP formation the nuclear modification 
factor at intermediate transverse momenta is dominated 
by the Cronin effect and at high transverse momenta by 
isospin effects and initial-state parton energy loss. Sur- 
prisingly, however, the contribution of final-state QGP ef- 
fects to the direct photon cross section in nucleus-nucleus 
collisions was small: less than —25% from the quench- 
ing of fragmentation photons, less than +25% from jet 
conversion and less than -f 10% from medium-induced 7 
bremsstrahlung. While experimental measurements are 
not yet precise enough to disentangle such modest ef- 
fects, they have already put severe constraints on theo- 
retical models that suggest a dominant role of jet-plasma 
interactions in the pt ^ 2 GeV part of phase space. In 
particular, PHENIX data is compatible with strong co- 
herent suppression of the bremsstrahlung 7, which is the 
main theoretical result of this paper. 

In conclusion, we suggest that only a systematic study 
of direct photons in heavy ion reactions for various sys- 
tem sizes and center of mass energies will help uncover 
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their full potential both as a baseline for jet tomogra- 
phy and as an independent probe of nuclear effects. In 
this exploration, precision d+A data is critical, since our 
theoretical results support the possibility that, in both 
proton- nucleus and nucleus-nucleus collisions, cold nu- 
clear matter effects play a dominant role in altering the 
cross section for direct photon production. 
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